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ABSTRACT 

We report the detection of rest-frame 6.2 and 7.7 /mi emission features arising 
from Polycyclic Aromatic Hydrocarbons (PAH) in the Spitzer/IRS spectrum of 
an infrared-luminous Lyman break galaxy at z=3.01. This is currently the high- 
est redshift galaxy where these PAH emission features have been detected. The 
total infrared luminosity inferred from the MIPS 24 fim and radio flux density 
is 2xl0 13 L , which qualifies this object as a so-called hyperluminous infrared 
galaxy (HyLIRG). However, unlike local HyLIRGs which are generally associated 
with QSO/AGNs and have weak or absent PAH emission features, this HyLIRG 
has very strong 6.2 and 7.7 /j,m PAH emission. We argue that intense star for- 
mation dominates the infrared emission of this source, although we cannot rule 
out the presence of a deeply obscured AGN. This LBG appears to be a distorted 
system in the HST ACS F606W and F814W images, possibly indicating that a 
significant merger or interaction is driving the large IR luminosity. 
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Subject headings: cosmology: observations — galaxies:dust emission — galaxies: 
high redshift — galaxies: mid-infrared 

1. Introduction 

Intensive star formation is a very common phenomenon in galaxies at high redshifts. 
Most high redshift galaxies are detected because of their high star formation rates, revealed 
by either the rest-frame UV luminosity in Lyman break galaxies (LBGs, Steidel et al. 2003) 
or at far-infrared wavelengths in, for example, submillimeter-selected galaxies (SMGs, Chap- 
man et al. 2003). It was long debated if LBGs and SMGs constituted different coeval pop- 
ulations. Recently Chapman et al. (2005) found that some SMGs have similar u — g and 
g — R colors as LBGs. The Spitzer Space Telescope permits exploration of the rest-frame 
IR properties of LBGs. For example, Reddy et al. (2006) studied the 24 /im properties of 
1.5 < z < 2.5 LBGs and considered their relation to other galaxies at that epoch, including 
SMGs. Huang et al. (2005) reported that 5% of the LBGs at z ~ 3 in the sample of Steidel et 
al. (2003) are detected at 24 /im and they defined this sub-population as Infrared Luminous 
Lyman Break Galaxies (ILLBG). ILLBGs have infrared properties very similar to sub-mm 
selected sources (Egami et al. 2004; Huang et al. 2005; Pope et al. 2006). Huang et al. used 
Spitzer [3.6]-[8.0] and [8.0]-[24] colors to argue that the 24 /im emission from the majority 
of ILLBGs is powered by intensive star formation. 

Observation of high redshift galaxies with the Infrared Spectrograph (IRS, Houck et al. 
2004) on board Spitzer permits the study of their 24 /im emission mechanism. Early IRS 
studies (Houck et al. 2005; Yan et al. 2005) have focused on 24 /im luminous (f24 >1 mJy) and 
optically faint sources. Most of the objects in both samples of Houck et al. (2005) and Yan 
et al. (2005) are identified as dusty AGNs at 1.5 < z < 2.5 whose spectra show power-law 
continua and strong silicate absorption but weak or absent emission features arising from the 
polycyclic aromatic hydrocarbon (PAH) molecules. Lutz et al. (2005) observed two lensed 
SCUBA sources at z ~ 2.8 with the IRS. Both sources are so-called Hyper-luminous IR 
galaxies (HyLIRGs, L 8 _i ooAtm > 10 13 L ) with strong PAH features. 

For galaxies at z > 2, the Spitzer MIPS 24 /im band detects rest-frame mid-IR contin- 
uum emission. At z — 2.1 and z = 2.8, the 6.2 and 7.7/im PAH emission features respectively 
enter into the 24 /im band. Thus, a high redshift, 24 /imselected-sample may preferentially 
select objects with strong emission-line features in the mid-IR bands. We have been studying 
just such a sample of LBGs, which have a redshift distribution centered at (z)=3.0 (Steidel 
et al. 2003). 
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This letter presents the IRS spectrum and multiwavelength properties of one ILLBG 
at z — 3.01 detected in the Extended Groth Strip (EGS) area. The multi-wavelength data 
set from All- wavelength Extended Groth Strip International Survey (AEGIS, Davis et al. 
2006) including ACS images permits the study of Spectral Energy Distribution (SED) and 
morphology of this source. We adopt Q m = 0.3, = 0.7, and H = 70 km s _1 Mpc -1 . 

2. Multi- Wavelength photometry and Infrared Spectroscopy of an ILLBG 

Most ILLBGs in the Huang et al. (2005) sample have 24 /im flux densities of ~0.1 mJy, 
too faint to permit efficient IRS observations. 1 Ashby et al. (2006) performed a wide-field 
deep u' and g' imaging survey with the MEGACAM on the MMT, covering the whole EGS 
and overlapping with deep Subaru R-band imaging and Spitzer/MIPS data. The goal of 
the MMT/Megacam survey is to obtain an ILLBG sample with higher 24 /im flux density 
for IRS spectroscopic follow-up. Huang et al. (2006) gave a detailed description of this 
sample. Objects with strong AGN activity such as QSOs are excluded from the sample 
using combinations of [3.6] — [8.0] and [8.0] — [24] colors (Ivison et al. 2004; Huang et al. 
2005). Spectroscopic redshifts for this sample are not yet available. We selected the object 
that is the subject of this letter, EGS20 J1418+5236 (hereafter as E21, the 21st in the IRS 
target list) as a LBG with u'-g'>1.4, g'-R=0.39, and R=24.30, and with very strong 24 /im 
emission, / 2 4 = 0.62 mJy. 

We obtained a spectrum for E21 with the IRS using the Low-resolution Long-wavelength 
(LL1) in staring mode, covering 20-39 /im. We observed the source for 24 cycles of 120 s 
duration. The total on-source integration time was 5851 s. We combined the reduced frames, 
subtracted residual background counts, and used the SMART package (Higdon et al. 2004) 
to extract calibrated one-dimensional spectra for both positive and negative beams. The 
final IRS spectrum is an average of both beams. The spectrum shows a significant detection 
of 6.2 and 7.7/xm PAH emission features at z — 3.01 ±0.016. This is consistent with a typical 
redshift for w-dropout-selected LBGs. Given the presence of multiple emission features in 
the IRS spectrum here, the redshift accuracy is significantly better than that reported for 
other z > 2 sources in the literature. Spectra of other objects often have low S/N, only a 
single emission feature, or featureless spectrum (e.g., Houck et al. 2005; Yan et al. 2005). 

The AEGIS multi- wavelength data set permits the study of E14 from X-ray to 20 cm 
radio band. While there is a single source in the MMT u', g', and Subaru R images, the 
HST/ACS images resolve this galaxy into two components, which we refer to as J 1 and J2 



see http://ssc.spitzer.caltech.edu/irs/documents/irs_ultradeep_memo.pdf 
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(Figure 2). In addition, another object, called J3, resides 3" to the east of Jl and J2, and 
could contaminate the Spitzer IR emission. Strong emission at u' from J3 suggests that it 
is at a much lower redshift, and using the full dataset we derive a photometric redshift of 
z = 0.79. J3 is not detected at 20 cm (5<x limit of /20cm <42 /iJy, compared to the 10cx 
detection of the J1+J2 component with /20cm =82 /iJy, Ivison et al. 2006). Therefore, 
we conclude that the J3 component does not contribute to the 24 /im flux density or the 
emission in the IRS spectrum. 

Both Jl and J2 are so-called u' dropouts (e.g., Steidel et al. 2003) and show low surface- 
brightness features that connect to Jl and J2 in the ACS V and I band images. The two 
components are separated by 0"8 or 6.24 kpc at z=3.01. These facts suggest that Jl and 
J2 belong to an interacting/merging system. The two components have different colors: Jl 
is red (Vab — I ab— 0.97), and J2 is blue (Vab — Iab—0.65). Jl is detected in both the J 
and K' bands, but J2 is not detected to Ki imity AB = 22 A. Jl has Jab — Kab=0.88 (or 1.9 
in the Vega system). We argue that Jl is mainly responsible for the IRAC, MIPS 24 /im, 
and radio flux densities and for the emission in the IRS spectrum. With Kab=20.89 for Jl 
and Kumit^AB = 22.4 for J2, J2 may contribute less than 11% of the infrared emission of this 
system (see Figure 3) . 



3. The Infrared Properties of E21 

Both the continuum and the 6.2 /im PAH emission feature contribute to the 24 /im 
flux density (/24=0.62 mJy) of E21. We measure the continuum emission as 0.55 mJy in 
the observed wavelength range of 20< A <32.5 /im, thus the contribution of the 6.2 /im 
PAH emission to the MIPS 24 /im flux density is 0.07 mJy. Most ILLBGs have 24 /im flux 
densities of ~0. 06-0.1 mJy (Huang et al. 2005). If other ILLBGs have comparable emission 
line fluxes, then this line would contribute significantly to their 24 /xm flux densities. 

The power mechanism for the mid-infrared emission of ILLBGs remains unclear. PAH 
emission features are associated with intensive star formation (Rigopoulou et al. 1999; Lu 
et al. 2003; Lutz et al. 2005). Rigopoulou et al. (1999) suggest that the 7.7/xm feature- 
to-continuum ratio, L/C, can be used to classify AGNs and starburst galaxies in ULIRG 
samples: L/C>1 for starburst galaxies and L/C<1 for AGNs. Though a few AGN dominated 
ULIRGs and most type-II AGN could have L/C>1 (Lu et al. 2003), a large 7.7/rni L/C value 
does at least suggest that intensive star formation is the dominant energy source. Peeters 
et al. (2004) argued that the 6.2/im PAH emission feature is a better indicator to separate 
starbursts and AGNs: QSO and type-I AGNs have typical feature-to-continuum ratio of 0.1, 
and starburst ULIRGs and type-II AGNs have a lower value in range of 0.5 < L/C < 2. The 
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6.2 and 7.7/zm feature-to-continuum ratios for E21 are 0.40 and 1.8 respectively, consistent 
with values of starburst ULIRGs. Furthermore, the limit on X-ray emission from E21 in the 
200ks Chandra observation (Nandra et al. 2006) is L 2 -ioKev < 1-0 x 10 43 ergs' 1 , suggesting 
the absence of strong AGN activity. Two other high redshift SMGs with IRS spectroscopy 
are identified as AGN dominated ULIRGs with similar IR luminosities and strong X-ray 
emission: SMM J02399-0136 at z=2.8 (Lutz et al. 2005), which much like this source has 
rather strong PAH emission; and CXO J141741. 9+582823 at z=1.15 (Le Floc'h et al. 2006), 
which has no PAH emission features. Furthermore, even when AGN are detected in high- 
redshift SMGs based on their X-ray emission, Alexander et al. (2005) argue that substantial 
ongoing star formation dominates the IR emission in these objects. Thus, we argue that 
intensive star formation is the dominant energy source for the infrared emission of E21. 
The different types of IR-luminous, high-redshift sources may represent different stages of 
HyLIRGs: starburst dominated E21 with strong PAH and no X-ray emission; intermediate 
type SMM J02399-0136 and other SMGs with both PAH and detected X-ray emission; AGN 
dominated CXO J141741. 9+582823 with no PAH but strong X-ray emission. 

Currently there is no sub-millimeter observation of this source. We use the 20cm radio 
flux and /sso^m/./^ocm-redshift relation (Carilli & Yun 2000; Dunne et al. 2000; Wang et al. 
2004a) to predict that E21 has an 850/im flux density of 10 < /sso^m < 17 mJy Similarity, we 
derive the rest-frame far-infrared luminosity using the FIR-radio luminosity relation (Condon 
1992; Carilli & Yun 2000): 

L FIR = 4nD^f 20cm W(l + z)-^ (1) 

-^8-lOOO^m — 

imL FIR (2) 

where q= log (/f/j?/ /20cm) =2. 34, and the radio spectral index a = —0.8. For /20cm — 0.082 
mJy at z=3.01, -^s-iooo^m is equal to 2xl0 13 L , qualifying this source as a HyLIRG. 

Fitting the SED with various galaxy templates yields a similar IR luminosity. We 
normalize the SEDs of two local ULIRGs, Arp 220 and Mrk 231, to the 24 fim flux density to 
predict the total infrared luminosity for this source: 6xl0 13 L Q and 2xl0 13 L respectively. 
(Fig. 3). The Arp 220 SED model predicts a much higher far-IR and radio flux density 
than observed. The Mrk231 SED provides a consistent interpretation to both the 24 //m 
and radio flux densities, as well as the 70 and 160 /xm flux limits. Similarly, figure 3 also 
compares the observed SED to models with IR luminosity L 8 _ 10 oo^m= 

10 i3.2-i3.3 Lg) (i.e., 

consistent with that derived from the radio emission) from Siebenmorgen & Kruegel (2006). 
The models broadly agree with the 24 /im flux density (and the limits at 70 and 160 /jm) 
and provide further evidence that the SED is dominated by emission from hot dust. This is 
also consistent with strong continuum emission in the IRS spectrum. Therefore, given the 
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constraints from the empirical and theoretical models, we conclude that hot dust dominates 
the far-IR emission in this object. 

This source is very unlikely to be amplified due to gravitational lensing. Gravitational 
lensing is very inefficient for J3 being a lensing galaxy at z=0.79. There are also no clusters 
or groups found in the area. Furthermore, the galaxy morphologies in the ACS image and 
color difference of Jl and J2 do not support strong lensing for this source. 

The morphology of E21 is reminiscent of those of local ULIRGs. Locally, about one third 
of ULIRGs have double nuclei with a mean linear separation of 6.2 kpc (e.g., Rigopoulou et 
al. 1999; Surace & Sanders 2000). E21 appears to have multiple components in the ACS 
images, see figure 2. The projected separation of Jl and J2 is 6.75 kpc. UV luminosity 
is a direct measurement for the relatively unobscured massive star formation in ULIRGs. 
Though ULIRGs are dusty, they typically have strong UV emission (e.g., Surace & Sanders 
2000), although this emission often stems from separate star-forming regions rather than 
those that dominate the IR (and bolometric) emission (e.g., Zhang et al. 2001; Wang et al. 
2004b). The U (A centra j=3414A) luminosities, vL v , for ULIRGs (Surace k Sanders 2000) 
cover a large range from 10 9 L Q to 5xlO n L Q with an average of 2xl0 10 L . At z~3, the 
observed J-band probes roughly the rest-frame U band. We use the J-band magnitude to 
estimate the rest-frame U luminosity for EGS20 J1418+5236, vL v {U) = 3.5 x 10 11 , roughly 
~l-2% of the bolometric emission. This is consistent with its local IR-luminous galxay 
counterparts. 



4. Conclusions 

IRS spectra of EGS20 J1418+5236, az = 3.01, u-dropout LBG in AEGIS show the 6.2 
and 7.7/zm PAH emission features. This is currently the highest redshift galaxy where the 
PAH emission features have been detected. This source has a total infrared luminosity of 
L 8 _ 10 oo/im=2x 10 13 L inferred from both the 24 /zm and radio flux densities. Nearby galaxies 
with these IR luminosities that are predominantly QSO/AGN, but given the presence of 
strong PAH emission features and the lack of any evidence for an AGN, we conclude that 
star-formation dominates the emission from this z=3 ILLBG. 

The 24 /zm-to-radio flux ratio, strong IRS continuum emission, and upper limits on the 
flux density at 70 and 160 /im implies that hot dust dominates the far-IR emission in this 
object, very similar to Mrk 231 and model SEDs of HyLIRGs. Interestingly, the theoretical 
models that best describe the shape and relative strength of the mid-IR SED require very 
high extinction values (A v > 70 mag). Because this object has strong UV emission, we 
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suggest that this LBG contains a deeply embedded starburst in addition to the relatively 
unobscured stars that dominate the rest-frame UV optical emission. This scenario may 
explain the general UV-IR properties of the ILLBG population. This example suggests that 
HyLIRGs at high redshifts can appear as either ILLBGs, or SMGs, or both. Larger samples 
of z ~ 3 galaxies with IRS spectroscopy and measurements of their far-IR SED are clearly 
needed to confirm these hypotheses. 

We wish to thank our collaborators on the AEGIS project for many insightful conversa- 
tions and much hard work. This work is based in part on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute 
of Technology under a contract with NASA. Support for this work was provided by NASA 
through an award issued by JPL/Caltech. 
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Fig. 1.— The IRS spectrum of EGS20 J1418+5236. The dashed line is the M82 SED 
shifted to z=3.01.The spectrum has remarkably similar 6.2 and 7.7 /zm PAH emission-feature 
strength and shape to those of M82. We cross-correlated the IRS spectrum with that of 
M82 to derive a redshift of z — 3.01 ± 0.016. The IRS spectrum at A >32/im becomes very 
noisy, so the 8.2/xm PAH emission feature is not detected. 
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Fig. 3. — The SEDs for both Jl and J2. The yellow boxes show the SED for Jl and 
cyan diamonds show the SED for J2. Upper limits on the flux-density (5er) are denoted 
by downward arrows. The blue solid line shows the SED for Arp 220 (Chary & Elbaz 
2001) and the red dashed line shows the SED of Mrk 231, in both cases normalized to the 
flux density at 24 /im. J2 is not detected in either the J- or K-band images. Because J2 
should not be redder than Jl, we use the Jl and J2 K-band flux ratio (with the limiting 
flux for J2) to estimate that J2 should contribute less than 11% to the IRAC and MIPS 
24 /im photometry. We use both the Arp220 and Mrk 231 SED models to estimate the total 
infrared luminosity of this source using the 24 /im flux density. The Mrk 231 SED model 
with L 8 _ 10 oo^m=2xl0 13 is in better agreement with the far-IR and Radio flux densities. 
Similarity, the gray-shaded region shows the range of theoretical SEDs of Siebenmorgen & 
Kruegel (2006) for starburst-powered IR-luminous galaxy models with Lg.iooo^m consistent 
with that estimated from the radio emission. Note that these models have not been rescaled 
to fit the observed photometry, but instead we show the predicted flux densities for a galaxy 
at z=3.01. The inset plot shows the IRS spectrum. The horizontal dashed line shows the 
observed 24 /im flux density (yellow box) and the other curves are the SEDs of Arp 220 and 

TVTrlr 931 sb for t.Vip main ficnirp 



